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Abstract 

We investigate numerically the quantum collision between a stable Helium nanodrop 
and an infinitely hard wall in one dimension. The scattering outcome is compared to 
the same event omitting the quantum pressure. Only the quantum process reflects the 
effect of diffraction of wave packets in space and time. 
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1 Introduction 



In recent years we have described the effect named: Diffraction of wave packets in space 
and time. The effect consists in a multiple peak wave train that is generated in a scattering 
event, by the interference between the incoming and the scattered packets. The wave train 
lasts to infinite time, due to the spreading of the incoming packet, that catches up with the 
scattered packet. [l]-[6] 

The effect occurs in matter wave scattering. It persists only for narrow enough packets 
compared to the scatterer extent. The effect was demonstrated numerically and analytically 
in various settings. As a possible laboratory experiment to verify the effect experimentally 
we suggested a collision between a Helium drop and a refiecting wall. [3] In ref.([3]) we 
omitted completely the self-interaction of the drop. The impinging packet that represented 
statistically the behavior of the drop was taken as a Gaussian packet, for which we could 
provide exact analytical expressions. However, the self-interaction provides the binding 
energy of the drop and hinders spreading, that is crucial for the appearance of the diffraction 
pattern. A more realistic calculation was needed. 

Helium is relatively a weakly bound liquid, even close to zero degrees Kelvin. Weaker 
than liquid water or other liquids at ambient temperature. The binding energy per Helium"^ 
atom in the superfluid phase is around Ei, — —7 .5^ K[19] whereas for liquid water it is around 
— SSOO'^A', depending on temperature and pressure. 

The Helium atoms are neutral and do not benefit from strong Hydrogen-like bonds for 
their binding. The apparent weakness of the interactions is nevertheless misleading. A 
Helium drop interacts with itself in a much stronger manner than a Bose-Einstein condensate 
(BEG) [7], and even there the nonhnear effects of the self-interaction are important. In 
particular the interatomic forces dominate over the so called 'quantum pressure' or quantum 
potential in the bulk of the drop. Inside the drop the density is constant and the quantum 
pressure vanishes. Near the edge of a drop, curvature effects take over and the quantum 
pressure becomes relevant. 

In the present work we investigate numerically the quantum scattering of a Helium droplet 
bouncing off an infinitely hard wall. The aim is twofold: Firstly an investigation like the 
present has never been carried out before and it is interesting in its own merit and, secondly 
we would like to have a more realistic treatment of the diffraction of wave packets effect in 
the context of liquid Helium. 

In order to simplify the treatment as much as possible we resort to the density functional 
approach. 

The density functional phenomenological method for interacting quantum systems, is 
based upon the solution of the one-body Schrddinger equation carrying self-interaction non- 
linear terms that depend on the density only. Many topics are treated nowadays by means 
of the density function theory. To name a few; electron transport in solids [8], electronic 
excitations [9], soft condensed matter[10] phase transitions in liquid crystals[ll], phonons 
in solids [12, 13], colloids[14], liquids and nuclei [15], atoms and molecules [16], quantum 
dots[17], etc. As a numerically viable phenomenological theory, is even becoming the dom- 
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inant method for treating many body quantum systems. This is evidenced by the number 
of papers on the topic that, by the late 1990's, surpassed the works using the Hartree-Fock 
method. [18] 

Density functional theories resort to the solution of a one-body nonlinear Schrodinger 
equation for a particle of mass m 

^^ = - + [Om') + U{v)]^l, (1) 

where 0(1^^1^) is a nonlinear and sometimes nonlocal[19] functional of the density p = 
and U{v) an external potential. The commonly used mean field equation for BEC systems 
belongs to this class of equations. [7] In this case it is referred to as the Gross-Pitaevskii 
equation. 

Density functional theory has been extremely successful in reproducing the properties 
of Helium nanodroplets.[19] Besides the mean binding energy per particle, average density, 
incompressibility and surface tension at zero temperature both for Helium^ and Helium^, 
in a stationary state[20], it accounts for capillary effects at low temperatures [21], the 
phase diagram of liquid- vapor coexistence [22], as well as the excitation spectrum (phonon- 
maxon-roton) of liquid Helium[19]. The model has been verified experimentally in scattering 
reactions [23]. It is currently being apphed to other areas of superfluid dynamics such as elec- 
tron bubbles, [24] adsorption on plates of Alkali metals[25], vortex line pinning[26], etc. 

The density functional theory is therefore a respectable method for the treatment of 
Helium nanodroplets scattering. In the next section we present results for the scattering of 
Helium nanodroplets in one dimension with and without the quantum potential. The last 
section provides some comments emerging from the observation of the numerical data. 
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2 A helium drop colliding with a wall 



A nonlinear, local, self-interaction of the Helium drop was proposed some time ago by 
Stringari and Treiner[20]. The parameters of the density functional were fitted to the known 
values of the binding energy per atom in the bulk, the incompressibility, the infinite atom 
matter density and the surface tension. Improvements to this functional seem to require 
nonlocal terms[7]. Such terms are engineered to reproduce the excitation spectrum of liquid 
Helium below the A transition point in a more accurate manner. In the present work we 
opt for the local version that is more manageable computationally and still quite successful 
phenomenologically. 

The operator 0{p) of eq.(l) is given by 



0{p) = 6p + ^(4 + 27)p^+'^ - 2 V'p (2) 

With h = -888.81 ^KX, c = 1.04554 10^ °K A^+^'f , 7 = 2.8, d = 2383°K 

For the onc-dimcnsional case treated presently we use the same parameter set, assuming 
independence from the y,z directions, i.e. an infinite planar slab. 

The stationary solutions of eq.(l) with the self-interaction of eq.(2), in one spatial dimen- 
sions arc obtained from the substitution ^(x,t) = e^^ '^(p{x), with p, the chemical potential. 
For finite size drops, p is a parameter , whose value determines the number of particles in 
the drop and tends to the bulk binding energy when the number of particles tends to infinity. 
For Helium"^, it is poo = —7.15^K} 

Stationary solutions are found by integration of the equation [20] 



. 2 

dp 



dpY 8rnp ( h p^ cp^. 

/^P+^ + ^^ (4) 



or equivalently 



dxj 8 mp d-\-l 

The integration of eq.(4), proceeds by a choice of p{x = 0) < p{x = 0)oo = 0.021836 A~^, 

the bulk helium liquid density for atomic matter. Eq.(4), with — |a,=o = 0, fixes the 

dx 

value of the chemical potential p > poo- Numerical roundoff error limits the possibilities 
of approaching p(x = 0)oo closer than 10^^ and still obtain a profile density that decays at 
infinite distances. We therefore chose the closest value we could and used this profile density 



^The transformation to standard energy units proceeds by multiplication with Boltzmann's constant k 
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Figure 1: Density profile of a Helium drop in A ^ of for p{x 
function of distance in A 



0) = 0.02183599 i'^ as a 



in all the scattering events to be presented below. The parameter for this density profile are 
p{x = 0) = 0.02183599, that yields p = — 6.71°/^. The number of particles per unit area 



N 



dx p{x) 



for this case is = 1.288 A ^, or an effective length X^fj 



N 



(5) 



59 A. The 



p{x = 0) 

drop then extends some 59 A along the x axis and is infinite along the y,z plane. The solution 
in the three dimensional case is not much different, in this case, the drop contains around 
4500 atoms. Figure 1 shows this initial density profile. 

The details of the scattering are not very sensitive to the number of particles and the 
extension of the drop. The majority of the particles reside within the bulk. The diffraction 
effects depend on the side wings of the distribution in figure 1 and are therefore almost the 
same for any number of particles, because the width of these wings is determined by the 
parameters of the nonlinear potential only. 

In the next section we will show results for smaller size drops whose profile resembles a 
Gaussian. For such drops the interference effect is cleaner, but not qualitatively different. 

The scattering starts with a drop impinging from the left onto an infinitely hard wall 
located to the far right. We took a drop centered at = 110 A with velocity v, defined by 



tp{x,t = 0) 



J, m V {x—xq) 



p{x) 



(6) 
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The wall is located at x = 150 A. We use a five step predictor corrector metliod[24] for the 

numerical integration. We demand an accuracy of more than 0.1% in both the conservation 

of normalization of the wave and the energy. This strict demand requires a tiny time step 

of around 10~^^ sec, and consequently lengthy runs of around 10^ iterations. A larger time 

step causes rapid deterioration of the accuracy and runaway behavior. A tiny time step is 

needed, due to the large values of the constants entering the self-interaction. 

Tn 

It was found that for packet velocities below around v = 50 — the drop collides with 

sec 

the wall elastically with a barely noticeable distortion. Such a threshold corresponds to a 
momentum transfer per particle at the wall of the order of Ap = 2 m v ^ 0.62A-\ This 
is a typical value for the excitation of ripplons, at the surface of a drop. [27] 

The drop is travelling without strictures, hence there is no superfiuid limiting velocity, 
it can remain superfiuid up to any velocity, in practice the drop fiows over a surface that 
limits the maximum frictionless speed. The experiment may also be viewed as one in which 
the wall is moved crashing into the stationary drop. We have performed such calculations in 
order to check our numerical schemes, and found them to agree with the results of a moving 
packet and a fixed wall as expected. 

Figure 2 exemplifies the result of the scattering for various velocities and times, such 
that Ax = V t, for all of them is the same in the free case. Above the threshold velocity 
there appears a multiple peak structure that receds faster than the packet. It resembles the 
diffraction in space and time structures seen in refs.[l]- [6]. For high enough velocities, the 
collision is inelastic. The energy is transferred into waves that recede faster than the bulk 
of the drop. The drop is almost stalled at the wall. In the bottom graph of figure 2 the 
collision process is still underway and the multiple peak structures are still being generated. 
There is a background hump under the peaks in the middle graph. This elevation is an 
incoherent background, similar to the one occurring in the diffraction in space and time 
phenomenon[l]-[6] for wide packets. In the next section we will show pictures for a thinner 
drop, for which the background is almost absent. 

Figure 3 shows the evolution of the scattering process for a velocity of v — 65.78 — 

sec 

The density fiuctuations inside the packet due to the collision, are progressively expelled 
out of the drop. This multiple peak tail wave, is a product of both the incoming and the 
refiected waves interfering with each other. 

In order to confirm the hypothesis that the interference effects are responsible for the 
multiple peak structures, we performed parallel calculations for the analog classical drop, by 
subtracting the quantum potential in the Schrodinger equation. In eq.(l) 0{p) is replaced 

by 0{p) - Uquantum, with Uquantum = " 2 m -^ "^ ' classical Scattering case , 

we use the same initial profile as the one depicted in figure 1, although it is not a stationary 
solution with the quantum potential subtracted. As a matter of fact there are no stationary 
solutions without the quantum potential. This may be easily seen by expanding around the 

^The awkward value of the velocity originates from a convenient choice of time units 
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Figure 3: Time snapshots of as a function of distance in units of A for v = 65.78— 
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Figure 4: as a function of distance in units of A for various velocities and times, classical 
case 

vacuum. 

Such an expansion produces the equation 



2 d 



P 



(7) 



whose solutions are oscillatory, with the density becoming eventually negative. 

Figure 4 depicts the classical scattering pictures corresponding to the quantum events of 
figure 2. As is evident from the figure, the classical collision does not produce a multiple 
peak tail. The collision results in large density fiuctuations in the bulk and noisy fluctuations 
at the edges. 

Tfl 

The corresponding classical packet evolution for v = 65.78 — appears in figure 5. 

sec 

Figure 6 shows classical and quantum scattering outcomes together at t = 60.8 psec. 
Clearly, the quantum scattering resists distortions in the bulk and gets rid of the energy in 
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Figure 5: Time snapshots of j^^l as a function of distance in units of A for v = 65.78;^, 
classical case 



10 





11 




12 



the form of coherent wave trains, whereas the classical scattering that lacks the quantum 
potential curvature-like term, renders the packet softer to deformations. The detail of the 
backwards region can be seen in figure 7. 

The coherence of the tail region is visualized in figure 8. In this picture we display the 
phase of the wave function for the scattering case of figure 6. The classical scattering phase at 
the bottom is coherent inside the drop corresponding to a travelling packet of fixed velocity, 
and incoherent or even random at the surface. The quantum phase is coherent both inside 
and outside (only part of the outside region is shown). The wavelength outside the bulk of 
the drop in the upper graph , for a; < 70 A is smaller than the wavelength in the bulk. The 
multiple peak tail recedes faster than the drop. 
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Figure 9: Density profile of a Helium drop in A ^ of for p{x = 0) = 0.02 A ^ as a function 
of distance in A 

3 Conclusions 

We have confirmed numerically that the effect of diffraction of wave packets in space and time 
does express itself in Helium^ nanodrop scattering. The results for a thick drop presented 
in the last section are not that clean as those of ref. [3]. A thick drop has a sizeable self- 
interaction and tends to conserve its shape in the collision. To enhance the effect, we 
considered also smaller size drops as the one depicted in figure 9. Here we used p(x = 
0) = 0.02 A~^ The number of particles per unit area is now = 0.26 and the 
effective extent is A^e// = 12.5 A. A three dimensional drop of this type contains around 
40 atoms. 

Figures 10 and 11 depict the classical and quantum scattering results for the profile of 

Tfl 

figure 9 and a velocity of f = 65.78 — at t = 60.8 psec. The asymptotic behavior of the 

sec 

quantum event is much cleaner than that of the thicker drop depicted in the previous section. 
In summary, it appears that it could be possible to observe the diffraction phenomenon 
described in [l]-[6] with Helium^ nanodrops by colliding them with a hard surface at high 
enough speed. 
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Figure 10: Quantum versus classical scattering for v = 65.78 — , p(x = 0) = 0.02 A 

sec 
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Figure 11: Quantum versus classical scattering for v = 65.78 — , p{x = 0) = 0.02 A 

sec 

enlargement of tail region. as a function of distance in units of A 
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